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Equilibrium thermal noise is known to destroy any quantum phase transition. What are the effects 
of non-equilibrium noise? In two recent papers we have considered the specific case of a resistively- 
shunted Josephson junction driven by 1// charge noise. At equilibrium, this system undergoes a 
sharp quantum phase transition at a critical value of the shunt resistance. By applying a real-time 
renormalization group (RG) approach, we found that the noise has three main effects: It shifts the 
phase transition, renormalizes the resistance, and generates an effective temperature. In this paper 
we explain how to understand these effects using simpler arguments, based on Kirchhoff laws and 
time-dependent perturbation theory. We also show how these effects modify physical observables 
and especially the current-voltage characteristic of the junction. In the appendix we describe two 
possible realizations of the model with ultracold atoms confined to one dimension. 
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I. INTRODUCTION 

Recent years have seen increasing interest in non- 
equilibrium many-body quantum systems. Schemati- 
cally, these systems can be divided in two categories: 
"closed", if isolated from the environment, and "open" if 
coupled to an external environment. In closed systems, a 
non-equilibrium situation can be generated by preparing 
the system in a given many-body quantum state and let- 
ting it evolve according to its (time-independent) Hamil- 
tonian. The resulting dynamics depends of course on 
both the initial state and the Hamiltonian, and is there- 
fore highly non-universal. To allow universal predictions, 
the initial state is often chosen as the ground state of 
some particular Hamiltonian. For instance, if this Hamil- 
tonian is close to a quantum critical point, the dynam- 
ics of the system is expected to follow universal scaling 
lawPEI. 

In open quantum systems, one is often not interested in 
the dynamics of the system, but rather in the properties 
of a non-equilibrium steady state, arising due to the flow 
of energy between two (or more) baths. In many cases, 
one bath is modeled in terms of a classical force, such as a 
voltage source, or an optical pump, while the second bath 
is treated in a full quantum manner. Depending on the 
details of the system, the quantum bath can be cither 
Markovian, as usually assumed in quantum optics, or 
non-Markovian, as more common in solid state devices. 
This difference gives rise to a different formalism needed: 
quantum master equation in the fo rme r case^Sl, and non- 
local real-time actions in the latteif^. 

Here we focus on non-Markovian quantum baths and, 
in particular, zero-temperature Gaussian baths. This 
type of baths can be obtained by integrating-out an infi- 
nite set of h armo nic oscillators, initially prepared in their 
ground stat d 9 * 10 *. In the absence of an external pump, this 
problem has been widely studied in the literature, espe- 
cially in the context of quantum phase transition a 11 " 15 !. 
The canonical example is a quantum particle in a double- 



well or in a periodic potential, under the effects of an 
Ohmic dissipative bath. If the coupling to the bath is 
weak, the particle occupies a coherent superposition of all 
minima of the potential. When the coupling to the bath 
becomes strong enough, the particle localizes in one min- 
imum of the potential, hence breaking the initial symme- 
try of the problem, through a universal quantum phase 
transition. In the case of a periodic potential, this transi- 
tion corresponds the insulator-superconductor quantum 
phase transition of a single resistively-shunted Josephson 
junction. The effects of an external drive on this systems 
is the subject of this work. 

In two recent papers^E2l ; we studied the steady state 
of a resistively-shunted Josephson junction driven by a 
stochastic voltage source, corresponding to 1// charge 
noise. Exploiting the scale invariance of the problem, 
we developed a novel analytic real-time renormalization 
group (RG) approach^HIID This approach offers a con- 
trolled way to describe the low-voltage properties of the 
junction. However, being expressed in terms of Keldysh 
path-integrals, it may appear highly non-transparent to 
the reader who is not familiar with this formalism. The 
goal of this paper is to study the same problem using 
simpler methods. Specifically, we will substantially relay 
on circuit theory, Kirchoff laws, and perturbation theory. 
Yet, we will reproduce all the main results obtained from 
the more involved RG calculations. 

This article is organized as follows. In Sec. |TT] 
we will introduce the specific model considered here, 
a resistively-shunted Josephson junction driven by 1// 
charge noise. The following three sections are devoted 
to the three main effects of the noise, namely: the shift 
of the transition (Sec. Ill), the renormalization of the 



resistance (Sec. |IV[), and the gene ration of an effective 
temperature (Sec. |V[). In Sec. VI we explain how these 
effects concur to determine the non linear current- voltage 



characteristic of the junction. Sec. VII concludes the ar- 
ticle with a brief summary and open questions. In the 
appendix we describe a possible realization of the model 
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FIG. 1: The noise-driven resistively shunted Josephson junc- 
tion: (a) Electric circuit of the junction, (b) Realizations with 
ultracold atoms (see Appendix) 

using ultracold atoms in one dimension. 



II. THE MODEL: A NOISE-DRIVEN 
SUPERCONDUCTING JUNCTION 

The object of the present study is the non-equilibrium 
device plotted in Fig. [T] As explained in the introduc- 
tion, the circuit consists of three main elements: 

1) A "pump". A stochastic, time-dependent volt- 
age source V^{t), capacitively coupled to the resistor. 
This voltage source models the so called "charge noise" , 
due to time-dependent fluctuations of charges in the 
substrate^. Being coupled linearly to the circuit, this 
type of noise has two main advantages: it can be easily 
introduced from the outside in a controlled experiment, 
and it allows for an exact analytical treatment as ex- 
plained below. In the following, we will consider in par- 
ticular the case of 1/f noise spectrum (Vn(uj)Vn{uj')) = 
i 7 o/(27r|ti;|)(5(a; — u/), where F measures the strength of 
the noise and has units of voltage-square. This parameter 
can be combined with the capacitance C and the electric 
charge (of a Cooper pair) 2e to form a unitless parameter 
F = F C 2 /(2e) 2 . As we will see this parameter plays an 
important role in physics of the problem. 

2) A "dissipative bath" . A linear resistor, correspond- 
ing to Ohmic dissipation, which is assumed to be at equi- 
librium (and specifically at zero temperature). The cou- 
pling to the bath is measured by the unitless ratio R/ Rq, 
where Rq — h/(2e) is the quantum resistance. The 
regimes R/Rq < 1 and R/Rq > 1 are often termed re- 
spectively "overdamped" and "underdamped" . The equi- 
librium quantum phase transition sits precisely at the 
boundary between these two regimes. 

3) A "quantum system" . A non-linear tunnel- 
ing junction, obeying the Josephson relations V(t) — 
(h/2e)dt4>(t) and I(t) — Jsin(0), where h is the Plank 
constant (in order to fulfill our didactic goal we shall 
not put h = 1 as often done), <p is the phase differ- 
ence across the junction and J is the Josephson cur- 
rent (proportional to the current of Cooper-pairs across 
the junction). The ratio between the Josephson energy 
Ej = hJ/(2e) and the charging energy E c = (2e) 2 /2C 
sets the third and last unitless parameter of the problem. 
The regime of Ej/Eq <C 1 is often termed "weak cou- 
pling" and is realized in ultra-small junctions (where the 



capacitance is strongly reduced, due to the small area), 
while Ej I Ec ^> 1 is termed "strong coupling" . 

The circuit of Fig[l] can also b e represented in terms 
of the quantum Hamiltoniar J 20 * 21 * 

H = + (2e)qV(t) - Ej cos(0) + H R [<f>) (1) 

Here the charge q is cononically conjugated to 4>. The 
last term Hr models the resistor as an infinite set of 
Harmonic oscillators. If the bath is initially prepared in 
the ground state, the corresponding degrees of freedom 
can be exactly integrated-out, leading to a real-time ac- 
tion with non-local kernels. To avoid this complication, 
here we will work directly with the original circuit, rather 
than with its equivalent Hamiltonian representation. 



III. RENORMALIZATION OF THE 
JOSEPHSON COUPLING 

As a first step, let us consider the weak coupling regime 
Ej <C E c and treat the Josephson coupling in a pertur- 
bative manner. If Ej = (J = 0) the fluctuations of 
the system are given by the linear sum of two terms (1) 
Johnson-Nyquist noise, due to the equilibrium fluctua- 
tions on the resistor; (2) external noise. In what follows 
we will assume that the two sources are statistically in- 
dependent. This is indeed a good approximation if the 
source of dissipation and the source of external noise are 
spatially separated (as in the case of Ref. |22j). 

The spectral properties of the equilibrium fluctuations 
of a resistor are well known and equal to (V^V*) = 
Rhujcotgh(huj/2T), where T is the temperature of the 
resistor. For later reference we note that, if T > 0, the 
low frequency tail of this spectrum is (V^V*) = 2RT, 
giving rise to white noise (delta-function correlated in 
time) commonly associated with classical thermal noise. 
In the zero temperature limit the spectrum becomes 

(KX) = Rh\w\ (2) 

This specific form of the spectrum, often called "quantum 
noise" , is a non- analytic function of u, giving rise to long- 
tailed correlations in time^. 

We now move to the voltage fluctuations induced by 
the external 1// noise. Perhaps surprisingly, we will 
find that the their spectrum has precisely the same fre- 
quency dependence as the quantum noise, giving rise to 
an interesting collaboration between classical and quan- 
tum noise sources. The origin of this behavior can be 
traced back to the RC circuit acting as a derivative of 
the incoming signalP^ More precisely, according to the 
linear circuit theory, the voltage over the resistor is — 
RI U = R/(R + i/ujC)Vn,uj- For frequencies significantly 
lower than 1/RC, one than obtains rj i{RC)ujVn ,u} ■ 
This approximation corresponds, in the renormalization 
group language, to the choice of an ultra-violet cutoff 
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(a) (b) (c) (d) 

R/R Q R/R Q R/R Q R/R Q 

FIG. 2: Schematic phase diagram of the noisy shunted Josephson junction: blue areas are superconducting and red areas are 
insulating, (a) Equilibrium phase diagram displaying a quantum phase transition at R = Rq. (b) First order corrections: 
the noise shifts the position of the transition in opposite directions, depending on the ratio between the Josephson coupling J 
and the ultraviolet cutoff Ao = 1/RC. The white area indicates the region where higher order corrections are expect to play 
a dominant role, (c) Second order corrections to the voltage response: the noise renormalizes the resistance of the junction 
and changes the slope of the transition line, (d) Second order corrections to the voltage fluctuations: the noise generates an 
effective temperature which transforms the transition into a universal crossover controlled by the non-equilibrium quantum 
critical point. 



Ao = 1/RC. For a 1// spectrum the voltage fluctua- 
tions: 

(V:V U ) = (RC)W(V^V N ^) = R 2 (2e) 2 ^-\u J \ (3) 

In the absence of the Jopsephson junction, the circui 
is linear, and the two sources of noise (equilibrium and 
non-equilibrium) simply sum-up: 

(V:V u )=Rh\L>\ + R 2 (2e) 2 ^\u>\ (4) 

Note that only the first component is multiplied by H, 
highlighting its quantum origin, while the second is of 
completely classical origin. 

A note of caution is now in place. Up to this point we 
have considered V as a classical field, giving rise to real 
correlations (V(t)V(t')}. To enable a quantum mechani- 
cal approach to the problem, one has to keep in mind that 
the expectation value S(t — t') = (V(t)V(t')) has both a 
real and imaginary part. At equilibrium, these two quan- 
tities are related by the fluctuation-dissipation theorem. 
In the presence of the noise, this relation is violated and 
the two quantities need to be determined independently. 
The real part S Rc (t-t') = i(V(t)V(t') + V(t'){t)) is re- 
lated to the fluctuations in the system and corresponds to 
a. The imaginary part, S Im {t - f ) = ±{[V(t), V{t')}), 
describes the response of the system to an external probe. 
For J = the system is linear and the response function 
is independent on the noise: Si m (u)) = cu^l/R + icuC)^ 1 . 

We are now in the position of adding back the Joseph- 
son coupling Jcos(0). For this task we need to compute 
the statistics of the phase fluctuations across the junc- 
tion. Using the Josephson relation V{t) — (h/2e)d t <j){t) 
and @ we obtain 

= 2 "^( 1 + i^ F ) log(1 + A ° (i_< ' )2) (5) 



Here we introduced back the capacitance, through the 
cutoff frequency Ao = 1/RC, in order to avoide the di- 
vergence of (pj at short times. 

Using eq.ip]) we can try to estimate whether the 
Josephson junction is capable of locking the phase across 
the junction <f> or not. Due to the weak dependence 
of the logarithm on its argument, we first estimate 
log(A |t-t'|) ~ 1 and obtain 5<j) 2 w R/Rq(1 + R/RqF). 
If the phase fluctuations are small, the Josephson cou- 
pling is able to localize the phase, driving the system 
towards a superconducting state. If, on the other hand, 
the phase fluctuations are large, the Josephson coupling 
will have nearly no effect. As a consequence, we may 
naively expect a transition at S(f> 2 — 1 or 

R* ( R* \ 
R Q V R Q J 

This handwaving argument can be substantiated 
through the powerful ideas of the renormalization group 
(RG) approach. In general, the goal of any RG method is 
to study the macroscopic behavior of a model by gradu- 
ally integrating over microscopic units. In our case there 
is no spatial dependance and the RG consists of averag- 
ing over fast processes. These fast processes renormalizc 
the tunneling coupling by "scrambling" the phase across 
the junction. To quantify this process, wc formally split 
the phase <j> into slow (s) and fast (f) components and 
average over the latter: 

Je** = Je^+'W -> Jfc-t'fye*** = J cS e^ (7) 

Here we used the property of Gaussian distributions, for 
which (e lA ) = e- 1 ' 2 ^. 

We now introduce an arbitrary frequency scale A sep- 
arating the fast from the slow processes. If we choose 
A J/(2e), it is reasonable to assume that these fast 
processes are independent on the Josephson coupling and 
their correlations are just given by Q. We than obtain 
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that 



6$ = (Mt)-MQ)T 



R A R r\ t j 1 

— 1 + — — F j du-r—r 

R Q \ R Q J J\u\>A M 



^( 1 + ^ F )'° g(i/a ° ) 

Combining Q and ^ one obtains: 



(8) 



JeS 

IT 



j 



A 
Ao 



-F -1 



(9) 



Eq.([9]) indicates the inset of a transition precisely at 
For R > R* the exponent in eq.([9]) is positive. 
Hence, as we reduce A towards zero the ratio J e g/A 
tends to zero. In this case, the Josephson coupling is 
said to be "irrelevant" and the junction behaves as an 
insulator. In the opposite case (R < R*) the exponent in 
|9]) is negative and the renormalized Josephson coupling 
grows, inducing a superconducting behavior. 

Until this point we considered only the weak cou- 
pling limit and performed a perturbative expansion in 
Ej/E c <C 1. We will now see that similar arguments ap- 
ply to the strong coupling limit Ej/E c 3> 1 as well. In 
this case, we can assume that the voltage across the junc- 
tion is zero. The current through the junction is given by 
two terms: quantum noise and external noise. In analogy 
to the previous analysis we find 



i? + F(2e)2 



Or, in terms of charge fluctuations, 



(Q(t)Q(t')) = (2e 



,2 R Q 

R 



R 
Rq 



F log 



t-t' 



(10) 



(11) 



From the above result we can guess the existance of a 
charge localization phase transition. If the charge fluc- 
tuations are larger than the charge of a single Cooper 
pair (2e) 2 , the junction behaves a superconductor. If the 
charge fluctuations are smaller, the capacitance becomes 
dominant and the junction behaves as an insulator. In- 
deed, an exact duality-transformatiorP^^ shows that the 
system undergoes a quantum phase transition at: 



R* 



1 



R* 
Rq 



F =1 



(12) 



Equation (|6| and ( 12 ) show the first non-trivial effect 
of the noise, namely the shift of the transition away from 
its universal equilibrium value R — Rq. This effect is 
schematically shown in Figj2^b). Note that in the weak 
coupling limit the critical resistance is smaller than Rq, 
indicating that the insulating regime is stabilized by the 
noise. On the other hand, in the strong coupling limit, 
the critical resistance R* is larger than Rq and the su- 
perconductor is stabilized. Our intuitive approach makes 
evident the origin of this effect: the noise increases the 



fluctuations of both the phase and the charge, always 
stabilizing the delocalized phase (i.e. the insulator at 
weak coupling and the superconductor at strong cou- 
pling). The same would be true for pure thermal fluc- 
tuations, which however would not change the critical 
correlations accordingly. 



IV. RESPONSE FUNCTION AND 
RENORMALIZATION OF THE RESISTANCE 

The second major effect discovered in Ref.pT] is the 
generation of an effective temperature. To understand 
the origin of this effect and its intuitive meaning, we 
need to consider higher order processes in the Joseph- 
son coupling, related to feed-back effects. As discussed 
above, the noisy RC circuit generates (non-equilibrium) 
voltage fluctuations with \uj\ spectrum, which translate 
into phase fluctuations, according to the Josephson volt- 
age law. Then, following the Josephson current law, 
the phase fluctuations translate into current fluctuations 
which are fed back into the RC circuit. When these cur- 
rent fluctuations pass through the resistor, they generates 
additional voltage fluctuations that correct the original 
voltage spectrum, and so on so forth (See Fig(3]) . 

To second order in the Josephson coupling we then 
have 



SS(t-t') = R 2 J 2 (sm(<j)(t))sm((f>(t'))) 



= R 2 J 2 e 



-i<(0(t)-0(t')) 2 > 



(13) 



As before, 5S(t — t') has both a real and an imaginary 
part, which renormalize respectively the fluctuations and 
the response of the system. To identify these two com- 
ponents it is useful to define two real functions, C{i) — 
Re((0(i) - <M0)) 2 ) and R(t) = lm((0(i) - </>(0)) 2 ) = 
^li[4>(t),(j)(t')]}, corresponding respectively to the corre- 
lations and response of </>(£)• For the specific choice of the 
cutoff introduced in RefP-2, C(t) is given by ( |5| a nd R(t) 
by R/ Rq atan(t/Ao). Using these definitions ~fl3[) we im- 
mediately obtain 5S(t) = SSn c (t) + i5Si m (t), where 



SS Rc {t) = R 2 J 2 cos{R{t))e- c{t) 
5Si m (t) = R 2 J 2 sin{R{t))e~ cw 



(14) 
(15) 



Let us first consider the contribution to the response of 
the junction SSi m (t). The Fourier transform of this func- 
tion is plotted in Fig. |4ja) at equilibrium (solid curve) 
and in the presence of a strong 1// noise (dashed curve). 
The comparison between the two curves shows a remark- 
able difference. At equilibrium the slope of the curve 
tends to zero as oj —¥ 0, while in the latter it tends to a 
constant. This trend is highlighted in FigQd), where the 
zero- frequency limit of the derivative of 5Si m (w) is shown 
as function of the noise strength F. At low frequencies 
we can approximate the contribution to the response of 
the system as 5Sj m ~ ojd u Si m (u! — 0). This term sums- 
up to the pre-existing (bare) response Si m w Ruj. We 



5 




Original voltage fluctuations: 
(Johnson-Nyquist and I/f noise) 



Renormalized 
voltage fluctuations 




V = RI 



Phase fluctuations 



Current fluctuations 



/=/sinW>) J K 



FIG. 3: Flow diagram of the feed-back loop generating a renormalization of the voltage fluctuations and the generation of an 
effective temperature (see text) 



conclude that the resistance is renormalized by 
dR{J) = d u S lm (io = 0) 



(16) 



The non-equilibrium renormalization of the resistance 
has important effects on the resulting phase diagram. 
Due to this term, the transition is now located at 
R + dR(J) = R* . Hence, with increasing J, the tran- 
sition moves to higher values of R (as shown in FigQb) , 
the correction is negative dR(J) < 0). When the Joseph- 
son coupling becomes of the order of the cutoff frequency 
(J/2e w A ) the transition is moved back to R ~ Rq. In 
the strong coupling regime, one should apply the dual- 
ity transformation described above. This analysis shows 
that, starting from J — > oo, when J is decreased the 
transition moves to lower values of R. The weak and 
strong coupling limit therefore lead to a consistent pic- 
ture, shown in Fig[2jc). 
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V. VOLTAGE FLUCTUATIONS AND 
EFFECTIVE TEMPERATURE 

We now move to the third main effect of the external 
noise, namely the generation of an effective temperature. 
For this task, we need to consider the contributions of 
the Josephson coupling to the voltage fluctuations of the 
junction S^ a (iS). The relevant expression is given in ( 14) 
and depicted in FigQc). We note a significant difference 
between the low-frequency behavior at equilibrium and in 
the presence of 1// noise. The equilibrium curve tends to 
zero at zero frequency, while the non-equilibrium curve 
has a finite zero frequency component. This feature is 
highlighted in Fig. Qb), where the value of the zero 
frequency component of 5Sn c is shown as function of the 
noise strength F. 

The zero frequency component of SSb, c is the origin of 
the finite effective temperature. As we discussed in the 
introduction, a finite temperature corresponds to white 
noise with spectrum 5r = RT. Here the spectrum is 
highly non linear, but at low enough frequencies, we can 
neglect the bare term ~ \uj\ and approximate S , r c (w) ~ 
5S-r,c(uj = 0) to obtain: 



R 



5Sr c (uj = 0) 



(17) 



FIG. 4: Second order contributions: (a) Renormalization of 
the voltage fluctuations SSb,c(<jj) at equilibrium (solid curve, 
R/Rq — 2, g/Ao = 0.1) and in the presence of 1// noise 
(red dashed curve, R/Rq = 2, F = 0.5, g/D = 0.1). (b) 
Zero frequency component as function of the noise strength 
F. (c) Renormalization of the response function <5i?i m (oj) at 
equilibrium (solid) and in the presence of 1// noise with F = 
0.5 (dashed), for R/Rq = 2. (d) Zero frequency derivative as 
function of the noise strength F. All axes are given in units 
of the cutoff frequency Aq = 1. 



It is worth noting that this effective temperature has to 
be understood only in an RG sense. It only determines 
the behavior of low-frequency behavior of the junction, 
while the high-frequency behavior is still strongly out of 
equilibrium. 

This effective temperature has drastic effects on the 
phase diagram, transforming the sharp phase transition 
into a smooth crossover. As we explained above, the 
first order predictions of a sharp transition were based 
on the power-law dependence of the renormalized Joseph- 
son coupling, eq. ([9]). However, as we now understand, 
this equation is valid only at frequency scales larger than 
T e ff/h. At this scale the thermal noise sets-in and leads 
to an exponential decay of the effective Josephson cou- 
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pling. For the superconducting behavior to be observable 
one has to require the renormalized Josephson coupling 
at the scale A = T e g to be larger than the effective tem- 
perature itself, or: 



J 



> 



(18) 



J 2 we find that, if R/Rq{1 



Recalling that T e g 
R/RqF) < 1/2, the condition ( 18 1 is always satisfied 
for J 0. In the intermediate regime 1/2 < R/Rq(l + 
R/RqF) < 1, on the other hand, a finite J is needed to 
obtain a superconductor. This effect is depicted schemat- 
ically in Figgd). 



VI. RENORMALIZATION OF THE 
CURRENT- VOLTAGE CHARACTERISTIC 

In the previous sections we described the effects of 
the noise on the renormalization of the model and on 
the resulting phase diagram. Here, we will describe how 
these effects can be probed by measuring the non-linear 
current-voltage characteristic of the model. Time depen- 
dent perturbation theory^ shows that the difference be- 
tween the total current in the junction and the current 
passing through the resistor I s = I — V/R is given by: 



I - 



V 
R 



= J 1 



J 

1 



dt ([cos ( <j>(t) 



2eVt 



,cos(<M0))]> 



h J 

dt e l(2e)yt// Tm(cos(0(i)) cos(0(O))) 



SSi m (u} = 2eV/h) 



(19) 



Here we used the definition of 6Si m given in eq.(15l 



Let us now study the behavior of this function in dif- 
ferent regimes. For small voltages, we can expand 5Si m 
in Taylor series of the frequency to obtain 



I = 



V 
R 



i - ^d u ss! m {u = o) 
.ft 



Inverting this expression we obtain: 
d u 5S lm (u, = 0) 111 " 1 



V/R 



sv 



= SR. 



(20) 



(21) 



This identity, already given in Sec |IV[ acquires now a 
clearer significance: the low-frequency slope of SSi m cor- 
responds to a renormalization of the low-voltage resis- 
tance. 

At larger voltages the junction deviates from this linear 
slope according to 



V 



1 



R + SR R 2 



[SW(u>) - LoduSimiuj)} 



(22) 



This quantity is plotted in Fig(5] on a log-log scale 
and clearly displaying an algebraic dependance. To 




FIG. 5: Non linear current- voltage characteristic, given by the 
difference between the supercurrent and the (renormalized) 
normal current I — V/(R + 5R), according to eq.(22l: equi- 
librium result (solid curve, R/Rq = 2) and non-equilibrium 
result (dashed curve, R/Rq = 2 and F = 0.5). The dotted 
lines represent the scaling result ~ y 2R / R Q-' 1 . 



understand this behavior, we consider the long-time 
limit of Sr m (*) - t -2R/R Q (i+FR/R Q )^ leading t0 

Sl m (t) ~ LJ 2R/R Q (l+FR/R Q )-l t Qr j _ V/(R + SR) - 

V 2R/k Q {x+FR/R Q )-x_ Thig power . law dependance is a 
direct consequence of the renormalization of the Joseph- 



son coupling discussed in Sec III 



Note that the non-linear behavior precisely disappears 
at the phase transition, where eq. (22) gives / ~ V. 



Thus, measuring the non-linear IV curve allows to deter- 
mine the precise position of the phase transition. This 
conclusion is however modified when one takes into ac- 
count the presence of a finite effective temperature. As 
discussed in Sec(Vj the predicted power-law behavior will 
terminate at the frequency scale u> = T ff. This poses a 
limitation on our capability of distinguishing between the 
different phases, thus transforming the sharp phase tran- 
sition into a smooth crossover. 



VII. SUMMARY AND DISCUSSION 

In this paper we consider the noise-driven resitively- 
shuntcd Josephson junction, originally proposed in 
Ref.[To] and [T7]. With respect to these two works, here 
we focus on the intuitive derivation and understanding 
of the results. Following the ideas proposed in Ref . [M] . 
we present the first order results of the renormalization 
group as a time-dependent average over fast processes. 
For the second order processes, involving the renormal- 
ization of the temperature and the resistance, we present 
even simpler calculations, based on time-dependent per- 
turbation theory. 

The resulting non-equilibrium effects significantly 
modify the resulting phase diagram, as probed by the 
current-voltage characteristic of the junction. The renor- 
malization of the resistance affects the slope at the low 
current limit of the curve. The renormalization of the 



FIG. 6: Effects of relevant and irrelevant noise sources, 
for R/Rq = 2. (i) Dashed curve: 1/f noise, (ii) Solid 
curve: a relevant noise source, entering the phase correlations 
through C(t) = ^log(l + A^ 2 +FA^ 4 ). (hi) Dashed- 
dotted curve: an irrelevant noise source entering the phase 
correlations through C(t) = log (l + A%t 2 + F-^\t\j . 



Josephson junction determines the non-linear behavior 
at larger voltages, and the effective temperature sets the 
transition frequency between these two regimes. 

One important effect which we did not consider here 
(nor in any of our previous papers) are the deviations 
from 1/f noise. It is known that experimental spectra 
always deviate from this theoretical curve. Based on di- 
mensional analysis it is natural to distinguish between 
noise sources affecting the long-time behavior of the cor- 
relations (relevant noise) and those which leave it un- 
changed (irrelevant). Without pretending to study this 
problem in depth, we have computed numerically the ef- 
fective temperature and the renormalization of the resis- 
tance for specific cases of relevant and irrelevant noise 
sources. 

In general, we found that the renormalization of the 
temperature is always present, independent of the rele- 
vance of the noise, in agreement with the general argu- 
ments given in SecjV] The renormalization of the resis- 
tance, on the other hand, strongly depends on the rele- 
vance of the noise. Two specific examples are given in 
Figj6] In the case of a relevant noise source (solid curve) 
the renormalization of the dissipation is extremely large 
even for pretty low noise strengths. In the case of an irrel- 
evant source (dotted-dashed curve), on the other hand, 
the renormalization of the dissipation is negligible. A 
complete understanding of these effects is still lacking. 
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for stimulating discussions. This research was supported 
in part by the US israel binational science foundation 
(EA and ED), the Israel Science foundation (EA), and 
the Swiss NSF under MaNEP and Division II (TG). 



FIG. 7: Experimental realization of the noisy shunted Joseph- 
son junction with ultracold atoms confined to one dimension. 
(See text.) 

Appendix A: Experimental realization with 
ultracold atoms 



In this appendix we discuss two possible ways to real- 
ize and probe the noisy shunted Josephson junction ([I]), 
using ultracold atoms confined to one dimension (See 
Fig{7]). In both realizations the dissipative bath corre- 
sponds to the low-energy excitations (phonons) of the 
liquid. In one dimension the phonons have universal 
properties^ their spectrum is linear and their spectral 
density constant. When integrated out, these modes pre- 
cisely correspond to a linear resistor^. 

The easiest way to understand the mapping between 
the Id system and the resistor is to compute the cor- 
relation function of the displacement field 6{x,t) in the 
Id model and to compare them with the (equilibrium) 
correlations of the phase across the junction 4>(t). If the 
phonons are prepared in their ground state, the correla- 
tions are: 



{6{x,t)6{x,t')) = ^(^)cos(^) 

f 7r/a 



(Al) 



K I dq — cos(cqt) « K log 

J-n/a M 



Here uj q ~ cq is the spectrum of the phonons and a the 
average atomic distance (which acts as a UV cutoff). In 
the second row we used the fact that the zero-point mo- 
tion of a harmonic oscillator is proportional to the inverse 
of its eigenfrequency and we introduced the proportion- 
ality constant K. For the proper choice of units, both the 
field 9 and the parameter K are unitless. We refer the 
reader to Ref|2H for a detailed derivation of the relation 
between K (the "Luttinger parameter") and the micro- 
scopic parameters of the models. At equilibrium (Al) 
and coincide, provided that we identify K — > R/Rq 
and 6(x, t) -> 4>(t). 

We now move to the physical realization of the Joseph- 
son coupling J. Following the pioneering work by Kane 
and Fisher^, our first proposed realization consists of 
a local impurity weakly coupled to the one dimensional 
system. As known from the literature^, the energy 
associated with a (back-scattering) impurity at posi- 
tion X Q is -ffimpurity = Vap(x = X ) = V cos(26> + 



8 



2ttXq/o). To drive the system out of equilibrium we 
propose to stochastically shift the position of the impu- 
rity as function of time, with 1/f spectrum, such that 
(Xo(lu)Xq(uj)) — (a/ir) 2 F/\u;\. To make a direct connec- 
tion with the noisy Josephson junction ([!]) it is enough 
to define a new variable <p(t) = \ (0(t) — irXo(t)/a). For 
this coordinate the bare fluctuations are given by the sum 
of an equilibrium (~ and a non-equilibrium (~ 

.F/|u;|) component and the non linear coupling (cos(0)) 
is time independent. 

Having established a formal equivalence between the 
noisy shunted Josephson junction and a vibrating im- 
purity in a one dimensional liquid, we now describe the 
physical consequences of this equivalence. For a shunted 
Josephson junction the natural physical quantity to look 
at is the non-linear FV curve. In the proposed realiza- 
tion, a finite current bias I = Q can be induced by drag- 
ging the impurity at a constant velocity (in addition to 
the random 1/f fluctuations), Xo(t) = al/(2e). The 
supercurrent I s ~ sin(</>) can be probed by measuring 
the atomic density at the distance a/2 from the impurity 
p(x = X +a/2) = cos(20 + 7rX o (t)/a + 7r/2), thus giving 
access to the non-linear I-V curve of the junction. 

Our second proposed realization is depicted in 
Fig{7^b). Here the model is realized in the anti- 



symmetric modes of two parallel Id liquids. If the tun- 
neling between the tubes is allowed only at a given po- 
sition, it immediately maps into a Josephson coupling 
^tunneling = t cos (5(f) , where 8<j) = 4>i — 4>2 is the phase 
difference between the two tubes and t/h the tunnel- 
ing rate. To introduce time-dependent noise in the sys- 
tem, we propose to apply a time-dependent potential- 
difference between the tubes, and to define accordingly 
(f> = 6(f> + J dt V(t). With respect to the previous realiza- 
tion, the voltage difference appears in the effective field 
with an additional integral over time. Thus, to mimic the 
1/f charge noise, we need to consider a voltage spectrum 
(V* (uj)V (lu)) ~ \w\. In fact, this type of noise is easier 
to generate than 1/f noise because its correlations decay 
algebraically rather than logarithmically. 

In this second realization, we can mimic a constant 
current bias by applying an additional DC potential dif- 
ference between the two tubes. The supercurrent can be 
probed by measuring the interference fringes of the two 
condensates (at the position of the tunneling junction), 
on the lines of Refl25l For high voltage differences the 
interference fringes completely disappear, corresponding 
to a linear I-V curve. As we lower the voltage, the fringes 
are expected to slowly reappear, indicating a non-linear 
I-V curve of the original model. 
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zo Apart from charge noise, superconducting junctions are 
affected also by "current noise" or "flux noise", due to 
time dependent fluctuations of the Josephson coupling. 
Fhis type of noise is highly non-linear, hard to predict, 
and will be neglected in this work. 

27 Recall that the derivatives of a function correspond to the 
momenta of its Fourier transform: if a function is non an- 
alytic, its Fourier transform must have long tails 

28 For instance, if the incoming voltage is a step function, 
the current through the resistor, and hence its voltage, has 
a short exponentially decaying pulse, similar to a delta 
function. 



